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Abstract 

We elaborate on the dynamics of ionized interstellar medium in the presence of 
hidden photon dark matter. Our main focus is the ultra-light regime, where the hidden 
photon mass is smaller than the plasma frequency in the Milky Way. We point out 
that as a result of the Galactic plasma shielding direct detection of ultra-light photons 
in this mass range is especially challenging. However, we demonstrate that ultra¬ 
light hidden photon dark matter provides a powerful heating source for the ionized 
interstellar medium. This results in a strong bound on the kinetic mixing between 
hidden and regular photons all the way down to the hidden photon masses of order 
lO " 20 eV. 



1 Introduction 


The nature of dark matter is one of the most exciting and pressing puzzles in particle 
physics and cosmology. Some dark matter candidates (most notably, the QCD axion and 
weakly interacting massive particles) have a more solid theoretical motivation than the oth¬ 
ers. However, given our current ignorance about composition and dynamics of the dark 
sector, a reasonable strategy in addressing the dark matter challenge is to leave no stone 
unturned, and explore as much consistent options and regions in the parameter space, as 
possible. 

One proposal in the off-the-beaten-path category, which attracted considerable attention 
recently (see, e.g., [1-6]), is that dark matter is made up of massive vector bosons, known as 
hidden photons. The simplest viable scenario for hidden photon dark matter assumes that 
no additional light particles carrying a hidden electric charge are present. In particular, this 
implies that the hidden photon mass term is of the Stuckelberg type, i.e. no Higgs field is 
involved. This is the setup we consider here. In addition to the hidden photon mass m, this 
minimal model is characterized by a single dimensionless parameter £, which determines the 
strength of kinetic mixing between the hidden and the ordinary photon [7]. It is convenient 
to work in the field basis such that the quadratic action is diagonal. The kinetic mixing 
translates then into a direct coupling between the hidden photon and the Standard Model 
(SM) electric current JF The resulting Lagrangian is 

£ = -4 - -F^ + —A^ - (1 + g2)1/2 ^ ( 4 . + eA.) > (!) 


where A p (A^) and F^ u (F IJU ) stand for the visible (hidden) photon gauge fields and field 
strengths respectively. Our notations are different from the standard ones, resulting in a 
somewhat unconventional factor of (1 + £ 2 )~ 1//2 in (1). Clearly, the two conventions agree 
at £ -C 1. The advantage of the one adopted here is that e is equal to the physical electric 
charge, as measured at high energies. 

The parameter space of this model has been constrained by a variety of laboratory exper¬ 
iments, astrophysical observations and cosmological arguments. In the present work we are 
interested in the ultra-light regime, m < ICC 11 eV. This region of parameters is especially 
interesting from the dynamical view-point because the interstellar medium cannot be treated 
as an empty space, as can be seen from comparing the hidden photon mass to the plasma 
frequency 
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7.4 • 10“ 12 eV 
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( 2 ) 


where n e ~ 0.04 cm -3 is a typical value of the free electron density in the vicinity of the Solar 
system. For many purposes plasma frequency may be thought of as an effective photon mass, 
so for m < uj p matter effects make the conventional photon heavier than the hidden one. If the 
hidden photon were heavier than the SM one in the intergalactic space, a resonant conversion 
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of hidden photons into SM photons would occur in the early Universe [2], Depending on 
the hidden photon mass, this may result either in cosmic microwave background (CMB) 
distortions or in the depletion of the amount of dark matter at low redshifts compared to its 
abundance at recombination. This results in the tight constraint on the ultra-light photon 
dark matter presented in Fig. 1. 

For m < 10" 14 eV the hidden photon is lighter than the SM one even in the intergalactic 
space, so the resonant condition is never satisfied and the bound disappears. The only 
existing constraint in this mass range comes from measuring the long-range magnetic field of 
Jupiter [8] and is not very restrictive (see Fig. 1). Existence of this constraint indicates that 
the plasma frequency is not identical to the photon mass. In particular, plasma frequency 
does not lead to screening of long-range static magnetic fields allowing to establish the bound. 

It is not a surprise that constraining ultra-light photons is hard. Indeed, in the limit of 
zero mass the Lagrangian (1) describes just the conventional electrodynamics with 

_ A u + eA^ 

M (1 + s 2 ) 1 / 2 

being the physical photon held, plus an additional massless vector held which is completely 
decoupled. In this case e is not a physical parameter and no constraint on its value is 
possible 1 . This demonstrates that extra care is required when discussing the ultra-light 
regime. In fact, additional rather tight constraints in the mass range 1CU 24 eV < m < 
10" 18 eV were claimed in literature (from the drift of fine structure constant induced through 
the Zeeman effect [1], and, most recently, from the Voyager magnetometric survey [9]). As 
we explain below, these constraints do not apply, because cosmic plasma shielding, neglected 
in [1,9], invalidates both effects. 

The principal goal of the present note is to improve on this situation and to push the 
bounds further into the ultra-light regime m <C 10" 14 eV. Not surprisingly, this requires 
a further careful look into plasma effects. The main idea is very simple. Hidden photon 
dark matter produces an oscillating electric field in the interstellar medium, which in turn 
induces electric current in its ionized component. Galactic plasma is not a perfect conductor, 
so the current dissipates resulting in a heat transfer from the hidden photon dark matter 
into ionized interstellar medium. By imposing that the rate of this heat transfer is smaller 
than the observed cooling rate from the dominant cooling mechanism for certain regions 
of the interstellar medium, we obtain constraints on the kinetic mixing coefficient e. The 
resulting constraint in presented in Fig. 1. It is somewhat stronger than the existing bound 
for m ~ 10" 11 eV. However, its most notable feature is that it extends deep into the ultra¬ 
light mass regime, all the way down to m 10" 20 eV. 

1 Of course, there is also no room for the hidden photon dark matter. 
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Figure 1: Bounds on the kinetic mixing £ as a function of the dark matter hidden 
photon mass m for three choices of parameters: black dashed lines corresponds to the 
cooling rate Q ~ 1CT 2 ' erg/s Uh, with nu ~ 0.1 cm -3 (an “optimistic” choice), red 
solid line to Q ~ 10~ 26 erg/s nn, with nn ~ 0.1 cm -3 (a “moderate” choice), and the 
blue dotted to Q ~ 10~ 26 erg/s nn, with nn ~ 1 cm -3 (a “conservative” choice) with 
T e = 8000 K, n e = 0.04 cm -3 in all the cases. The two filled regions show the CMB 
bound [2] and the bound from the Jupiter magnetic held [8]. This plot does not include 
constraints from black hole superradiance, which are briefly discussed in Section 5. 

2 Plasma equations with hidden photon 

To derive this result let us study the dynamics of the photon and hidden photon fields in 
the interstellar medium, which we model as a non-relativistic plasma. Then the motion of 
electrons and ions is described by the momentum conservation equations (see, e.g., [10]) 


du e _ 

m^ + meK-VK 


m ,T- - m . ( v , ■ V | a . 


e 

+ sE + u e 
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where tq and u e are the velocity fields of the electrons and ions respectively. The first term 
on the r.li.s. of both equations represent the Lorentz forces experienced by charged particles 
due to the conventional electromagnetic fields and due to the hidden one. The last term 
on the r.li.s. of both equations accounts for the friction due to electron-ion collisions. The 
frequency of electron-ion collisions v is given by 


v = 


Ay/2na 2 n e n / n e \ /8000 K\ 3//2 

^w^ logAcK3 ' 10 ' eV (ao4==0 —) 


( 5 ) 


where n e is the electron density, T e is the electron temperature and the Coulomb logarithm 
is given by 

/ 47tT 3 \ 1/2 

log Ac = log ( J ~25, 

V OL 6 n e J 


for typical parameters of interstellar medium (n e ~ 10 2 cm 3 , T e ~ 8000 K). 
Eqs. (3) are complemented by the electromagnetic and Proca field equations 


HA* = 


DA 11 + m 2 A» = 


(1 + e 2 ) 1 / 2 

ee 


J* 


(1 + e 2 ) 1 / 2 


J* 


and the fields satisfy 

d,,^ 1 = d.A 11 = 0 , 

as a consequence of the Landau gauge choice for the electromagnetic field and the Proca con¬ 
straint for the hidden photon field. Finally the current is written in terms of the electron/ion 
density and velocity fields as 


J» = 



(rq, UjMj) (n e , n e u e ) 


( 6 ) 


To study wave propagation in the cold plasma we linearize the above system of equations 
around a configuration with zero photon and hidden photon fields, zero electron velocity and 
constant electron and ion densities n e and rq. That is, we linearize the system by considering 
the (constant) densities to be zero-order quantities and the velocity fields and electric and 
hidden electric fields to be first-order quantities. The magnetic fields can be neglected at 
the first order. The linearized momentum conservation equations take the form 


m e 


du e 

~di 



c -> 

- (1 + £ 2y/z ( E + eE ^> - meV e _ “0 

6 z* 

(1 +£ 2y/2 ^ E + SE ) _ ~ U e ) ■ 


( 7 ) 

( 8 ) 


At this order the total momentum of the system is conserved at each point, so that in the 
rest frame we find 

m e u e + iriiUi = 0 
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at all times. Then at the leading order in m e /rrii we may set the ion velocity to zero, Eq. (8) 
holds trivially and we are left with Eq. (7) at v n — 0. 


3 Plasma shielding 


After switching to the Fourier modes for all first-order quantities, 

A(x,t) = 1 /2 J d 3 kdcuA(k,uj)e isn - iut (9) 

we may solve from Eq.(7) for the electron velocity in terms of the electric fields as 


u e (k,uo) = 


E + eE 


m e {y — ioj) (1 + £ 2 )V 2 ' 


( 10 ) 


Plugging this into the electromagnetic equations for the photon and hidden photon fields 
and decomposing the spatial parts of the vector fields into transverse and longitudinal parts 


Ai(k,u>) = —b Ati i Axiki = 0 

k 2 


we find 


k 2 -u 2 + ( 1+e 2 ) 


2\-l 


Q2 

“p 




H 2 £ + m 2 (l + e 2 ) 
for the transverse components of the vector fields, and 


n 2 
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a 2 p s 


1-i 


-U + (1 + £ ) I n 2 £ n 2 £ 2 + m*(l +A) 
for the longitudinal component, where 


r>2 _ 02 p 

p ~ 1+* 


A 1 

Ai 


= 0 


A, 


= 0 


(11) 

( 12 ) 


(13) 


We see that the presence of matter affects the vacuum mixing of the two fields A and A, 
as anticipated. We may now diagonalize the previous set of equations to obtain the two prop¬ 
agating modes and their respective dispersion relations for the transverse and longitudinal 
parts of the fields. 

Let us first discuss the structure of the propagating modes. For the purpose of the dark 
matter discussion we are interested in the non-relativistic regime k oj. Then the transverse 
and longitudinal equations approach the same form and the propagating modes correspond 
to the eigen-modes of the matrix 


M 2 


(1 + £ 2 ) -1 


n 2 n 2 e 

n 2 e n 2 p £ 2 + m 2 ( i + £ 2 ) 


(14) 
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At 77i <C Up the two eigenvalues of the matrix M 2 take the form 


cof, = ni[i + o[ — 


m 


771 


UJu = 


l + £ 2 


. 771 

1 + 01 (P 
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(15) 

(16) 


Ultra-light dark matter is comprised of the second (lighter) mode A^m- The corresponding 
A and A components take the form 


A — —Adm 
A = A dm 


( 1 + E 2 ) 1 / 2 

1 


1 + 


m 


fi 2 (l + £ 2 ) 


m 2 £ 2 


(1 + £ 2 )V2^ H 2 (l + £ 2 ) 




( 17 ) 

(18) 


To see the plasma shielding effect let us calculate the corresponding observable electromag¬ 
netic held given by 


A 0 bs 


A T eA 


e m 


A 


DM- 


(19) 


(1+e 2 ) 1 / 2 l+e 2 !J 2 

We see that on top of a vacuum mixing suppression present at small £, interactions of ultra¬ 
light hidden photon dark matter with the SM matter are suppressed by an additional factor 


FT/ % for i/ « m « Wj, 

O 2 ~ I 051 (,„■ „ ' 

It is straightforward to understand the physical origin of this suppression. Given that the 
plasma frequency is much higher than the hidden photon mass the light mode should be 
highly aligned with a linear combination of A and A, which does not interact with plasma, 
in order to stay light. As a consequence of this suppression, direct detection of ultra-light 
photons becomes very challenging. For instance, strong bounds on the value of £ in the 
mass range 1CU 23 eV< m < 10 -18 eV were claimed recently in [9], based on the precision 
magnetometry performed on board of the Voyager space probes. However, including the 
plasma shielding effect derived here and neglected in [9] , reduces the expected signal at least 
by a factor of 10 -12 (for m = 10“ 18 eV and taking n e ~ 10“ 3 cm -3 , which is appropriate for 
the heliosheath region [11]) completely eliminating the constraint. 

The plasma shielding suppression cannot be directly applied to the Earth based experi¬ 
ments, in particular to the constraint based on the Zeeman effect [1], and to the SQUID based 
setups proposed in [9], because these experiments are performed in the dielectric environ¬ 
ment of the Earth atmosphere. However, a strong signal suppression of a very similar origin 
applies in this case as well. Indeed, here we encounter the situation analyzed in detail in [5]. 
Namely, effectively the experiment is performed in a vacuum (dielectric) cavity surrounded 


6 









by a conducting medium. As demonstrated in [5], for small photon masses the electric field 
inside the cavity is suppressed then by an additional factor of ma x{0(m 2 R 2 ),0(mRvDM)} 
and the magnetic field is suppressed by an additional factor of 0(mR), where vjjm ~ 1CT 3 
is the dark matter velocity and R is the size of a cavity. For the mass range m < 1CD 19 eV 
considered in [1] the bound disappears even if one estimates the size of a cavity to be 
R ~ REarth ~ 6000 km. Given that the Earth soil has very decent conducting properties it 
is likely that the actual suppression is much stronger and corresponds to R ~ 60 km, which 
is the height of the Earth ionosphere. A proper accounting for this suppression is definitely 
required to evaluate the feasibility of the SQUID magnetometry proposal of [9]. 

To conclude this discussion it is worth noting that the plasma shielding effect does not 
affect the bound based on the study of the static magnetic field of the Jupiter. Indeed, this 
bound is based on the study of a kinematical regime k oj, which is very different from the 
non-relativistic regime discussed so far. Plasma does not screen static magnetic fields and 
the Jupiter bound applies. 


4 Heating up the Galaxy 


To summarize the above discussion, the plasma shielding effect leaves hidden photon dark 
matter practically unconstrained for the mass range 2 1CU 23 eV < m < 10~ 14 eV. As we will 
see now, a further look into plasma effects allows nevertheless to establish a quite stringent 
constraint over a sizable part of this mass range. The idea behind this constraint is that 
an oscillating electric held induced by hidden photon dark matter should dissipate as a 
consequence of non-zero resistivity of the interstellar plasma. The resulting heat transfer 
rate from dark matter into the interstellar medium should not be too fast. 

To calculate the heat transfer rate we need to End the imaginary part of the frequency 
corresponding to the dark matter mode. As before, we restrict to the non-relativistic limit, 
so that the eigenfrequency equations are the same for the transverse and longitudinal modes. 
Even though we are mostly interested in the ultra-light regime, m -C oo p , for completeness 
we consider also the opposite regime, m uj p . The dispersion relations of non-relativistic 
plasma waves are determined from the following equation, resulting from (12) (or, equiva¬ 
lently, from (13)), 


u = 



cot, 


(i + v) ± \ 


rrr + 




4 ujp r m 2 


-\ 


(l + i£) J (l+-V>(l + t) 


( 21 ) 


2 Barring the Jupiter constraint for the heavy end of this mass interval, which is rather mild. See also 

section 5 for a brief discussion of the superradiance constraint. For masses much lighter than 10~ 23 eV, the 
hidden photon is not a viable dark matter candidate, because its de Broigle wavelength is longer than the 
observed small scale structures. 
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As before, the two branches correspond to the conventional plasma oscillations and to the 
dark matter mode. They exchange their roles for m above and below c o p . Namely, for m <C c <j p 
it is the negative branch that behaves like the hidden photon, whereas for m 3> co p it is the 
other way around. Taking this into account and decomposing the frequency into its real and 
imaginary parts oj = Uh + i^h we have in the two limits 


f yyyr for m < u p 
\ m 2 for m cu p 


( 22 ) 


and 


for rn < u p 

for m » Up 


(23) 


From this last expression we obtain the rate of the energy loss for the dark matter hidden 
photon field 

Q = 2|q| Ph , 


where 


in the 

Q ~ 


Ph is the dark photon energy density. This translates into 

n ~ 2 • 10-V^ ( _ - _V ( _^_ -) 

scm 3 \7.4 • 10 -12 eV/ v0.3 GeVcm -3 / 
light mass regime (m -C c o p ) and into 

2 ■ ( 7 - 4 ' lirl2eV V ( p ; , ) (^ 

scm 4 \ m J VO.3 GeV cm 6 J \ l e 


^8000 K^ 3/2 



n P 


0.04 cm -3 


2 


(24) 


(25) 


in the heavy mass regime (m uj p ), where in both cases we assumed e 2 1. Interestingly, 
in the light mass regime the energy loss rate is independent of the electron density. This 
comes out as a result of the cancelation between the two competing effects—decreasing the 
electron density suppresses the collision rate z/, but at the same time the plasma shielding 
effect is less pronounced in the dilute plasma. 

This energy gets dumped into a thermal motion of plasma electrons. The current state 
of the interstellar medium is maintained through a delicate balance between various cooling 
and heating mechanisms (see, e.g., [12,13] for an introduction), which will be destabilized if 
the heat transfer rate from the dark matter photons is too fast. Unfortunately, there is no 
direct measurement of the heating rate, although theoretical expectations for a number of 
key processes are available. On the other hand, there are more observational handles on the 
cooling rate (and it is encouraging that the results match well theoretical expectations for 
the heating rates). So to obtain a bound on the hidden photon parameters, we impose that 
the rate given by (24), (25) is smaller than the observed cooling rate for various regions of 
the interstellar medium. 

One way to arrive at the bound is to consider the warm ionized medium, a dilute phase 
of the interstellar medium of temperature around 8000 K [13], extending about the galactic 














midplane with a thickness of more than 1 kpc, filling a volume fraction / 0.4 ~7~ 0.2 [14]- 

For this phase, the cooling is dominated by the inelastic collisions of electrons and hydrogen 
atoms with singly ionized carbon atoms Cn [13]. The ground state of Cjj is split into two 
fine structure levels, and inelastic collisions can excite transitions from the 2 P \/2 level to the 
2 F 3 /2 level (excited atoms being labeled C* n ). Subsequent de-excitation of the excited fine 
structure level yields the far-infrared emission in the 157.7 p,m line. 

The rate of cooling due to this mechanism can be deduced from the abundance of C* n in 
the warm ionized medium. For different regions of the interstellar medium this abundance 
can be approximated from the integrated column density of C* n along different lines of sight. 
The latter can be in turn obtained from the C* n absorption lines in the far ultraviolet given 
by transitions originating in the 2 P 3 /2 level. 

The procedure outlined above has been carried out by several authors [15-17] to obtain 
the cooling rate from Cjj emission into the infrared. In [17], measurements of the 1037.018 A 
and 1335.708 A far UV absorption lines of 43 objects at high galactic latitude (b > 30°) are 
used to obtain the cooling rates per hydrogen atom for the clouds along the corresponding 
lines of sight. The gas probed along these lines of sight is claimed to be at least partialy 
in the warm ionized phase. The results for the cooling rates vary between an average of 
~ 10” 26 n# ergs/s for directions corresponding to low velocity clouds and a smaller value of 
around ~ 10 ~ 27 nn ergs/s for higher velocity clouds, where n# is the density of Hydrogen 
atoms. The latter is itself expected to be in the 0.1 -t-1 cm -3 range [18,19] (and may also be 
estimated from the data presented in [17] from the available values of the electron density 
and ionization fraction). Imposing the condition that the heating rates (24) and (25) are 
smaller than these cooling rates we find the bound presented in Fig. 1. We accounted for the 
uncertainty in the value of rifj and other parameters, by presenting the bound for several 
choices of the parameters, as explained in the caption. 

5 Conclusions 

To summarize, in this note we elaborated on plasma effects in the ultra-light hidden photon 
dark matter model. On one side, we pointed out that the previously neglected plasma 
shielding effect makes it quite challenging to directly detect ultra-light hidden photon dark 
matter. On the other hand we derived a quite stringent bound on the mixing parameter 
of ultra-light hidden photon dark matter, resulting from the measured cooling rates in the 
Galactic warm ionized medium. It appears plausible that our bound may be somewhat 
sharpened by a dedicated analysis of the cooling rate in specific regions or in a different 
environment. Given that the heat transfer rate from dark matter to the interstellar medium 
is independent of the electron density, one may consider, for example, also warm HI regions. 
A brief look at the corresponding typical cooling rates and ionization fractions presented 
in [12] suggests that indeed this results in a comparable bound. It would be interesting to 
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see whether the bound can be further improved by considering the intergalactic medium 
either in the current epoch or at higher redshifts. The heat transfer rate in this case gets 
suppressed compared to what one finds in the Galaxy due to smaller dark matter density. 
However, this may well be compensated by an even stronger decrease in the cooling rate, 
which is proportional to n 2 . We leave this kind of analysis, requiring more detailed study of 
astronomical data, for the future. 

Our bound becomes weaker as the photon mass decreases, and disappears at masses 
around 10~ 2 ° eV. It will be interesting to explore other ways to constrain ultra-light hidden 
photon dark matter with these tiny masses. We are aware of two ideas how this can be 
achieved. First, similarly to ultra-light axions [20,21], also ultra-light vectors in the mass 
range ~ 10~ 22 -3 10~ 10 eV may be constrained by measuring spins of astrophysical black 
holes, through the effect of black hole superradiance. Existing spin measurements for stellar 
mass black holes constrain the mass range rsj 10 12 -T -10 11 eV [22], It is unclear whether the 
present spin determinations for supermassive black holes are robust enough to use them for 
establishing the bound. By taking them at the face value, one may exclude the mass range 
~ 10 -20 -310~ 16 eV [22,23]. Note, that the superradiance bound does not rely on a particle 
being a significant part of dark matter, and applies at e = 0. On the other hand, matter 
effects may inhibit the superradiance instability at sufficiently large values of e. To the best 
of our knowledge, these were not studied for hidden photons, so at the moment it is unclear 
what is the largest value of £ such that the bound applies. 

In addition, it was proposed recently in [24] that ultra-light scalar dark matter in the 
mass range ICO 23 -310~ 22 eV produces an oscillating gravitational potential, which may be 
observed with pulsar timing arrays. It is natural to expect that this also applies to the 
ultra-light photons in the same mass range. It will be interesting to study whether the two 
models can be distinguished observationally. 

We would like to thank Andrei Gruzinov for very helpful discussions. This work was 
supported in part by the NSF CAREER award PHY-1352119. 
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